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A motion of neutral and charged particles located nearby a metal surface under joint 
action of time-periodic electric field E(t) directed normally to the surface and perma- 
nent magnetic field H directed along the surface has been considered. It has been shown 
that due to a nonlinearity of the system there exists a directed transport of the particles 
perpendicularly to both magnetic and electric fields with velocity which is proportional 
to E 2 H. The velocity has been evaluated for charged (ions, electrons etc.) and neutral 
(atoms, molecules, nano-size clusters etc.) particles adsorbed on a metal surface. Cor- 
responding surface electric current has been found as a function of a frequency of the 
electric field and material parameters. It has been noted that the similar phenomena 
appear in a bulk of electrically inhomogeneous media. 



PACS: 05.45.-a, 05.60.-k 

The directed transport of particles under ac- 
tion of a time-periodic force is now well es- 
tablished phenomenon in various non-linear sys- 
tems. The details of this transport has been stud- 
ied for one-dimensional motion of particle in the 
space-periodic potential U(x) and a time peri- 
odic force having general form F(t) = Fo costut + 
Fi cos (2ut + 9). (See review and recent works 

H . ) Main physical reason for such directed mo- 
tion of the particles is due to the so called la- 
tent asymmetry of driving force F(t). Though the 
time average of F(t) is zero the time average F 3 (t) 
doesn't vanish and depends on 9. Nonlinearity of 
the system leads to mixing of harmonics such a 
way that average velocity is not zero and directed 
along F 3 (t). 

The aim of this paper is to show that the similar 
time- asymmetry arises in a permanent magnetic 
field even in the case of purely harmonic external 
force F(t). The phenomenon is closely related with 
a conventional Hall effect in conducting media. In 
difference with the latter the DC electric current is 
arising on interface under a joint action of the har- 
monic electric field directed normally to the surface 



and a permanent magnetic field directed along the 
surface. The effect is of the second order with re- 
spect to the amplitude of the electric field. For this 
reason we'll call it the nonlinear Hall (NHL) effect. 

Firstly we calculate the transverse velocity of a 
single charged particle located on a metallic surface 
in the AC electric field and longitudinal permanent 
magnetic field. Then, in Sections 2,3 we consider 
the similar phenomenon in a surface electron gas of 
metals and semiconductors. In Section 4 we eval- 
uate the transport velocity of neutral particles in 
above conditions. Finally we give a phenomeno- 
logical description of the NLH effect in inhomoge- 
neous conducting media (Section 5). 

I. DIRECTED MOTION OF SINGLE 
CHARGED PARTICLE. 

In order to demonstrate the physical mechanism 
of DC current let us consider a motion of a single 
ion with a mass M on metallic surface. A complete 
description of an interaction between the ion and 
metallic surface is a quite complicated problem. In 
framework of a density functional approach it has 
been discussed in a recent monograph by Liebsch 
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H . There are a few solidly established facts we are 
going to use in our paper. The interaction between 
the ion and the metallic surface consists mainly of 
two parts. The first part is an attractive Coulomb 
force between the ion and an image charge having 
a potential Vi m (z) 



V im (z) = -Q 2 /(2ze) 



(1) 



where z is a distance between the ion and an edge 
of the metal surface and e is a static dielectric con- 
stant. We gave here the simplest expression for 
such force at large z neglecting the dynamics of 
the metal electrons. More rigorous consideration 
includes an integration over the complex surface 
polarisability, which modifies the z- dependence of 
Vi m {z). For the sake of simplicity we will use (1) 
with some dimensionless parameter c\ 

V lm {z) = - Cl Q 2 /z (2) 

The second part of the interaction is connected 
with a short-range repulsive exchange potential 
V ex (z) which we choose in exponential form 



V ex {z) = c 2 aiexp{-(> 



(3) 



were z e is an equilibrium distance between the phy- 
sisorbed ion and metal surface, a\ is a radius of 
exchange forces on surface. The c 2 is chosen from 
a equilibrium condition 



V{z) = V im {z) + V eso {z) 



dV/dz\ 



c 2 = cxQ 2 /zl 



(4) 



(5) 



(6) 



The motion of the physisorbed ion with a velocity 
v leads to electronic excitations in the metal and 
therefore to a friction force for such motion. 
The friction force may be written as |q] 



Here v is the particle velocity and m, r/± are the 
friction coefficients for the motion parallel and nor- 
mal to the surface. If the charge is located outside 
the electronic density of metal, it can be shown 
that 77|| = 77j_/2 ||. At large ion surface distances 
z the friction parameter v± depends on z as 

Q 2 uj f 1 



M k F w 2 l z 3 



(8) 



where hdjp,kp,l and lo p are a Fermi energy, a 
Fermi momentum, a mean free path of electrons 
and a plasmon frequency respectively. 

For rough estimation it is convenient to use more 
simple formula 



Q 2 



(9) 



-M{ri\\i>\\ + j]±v± 



(7) 



Ma z 3 

Here <7o is a bulk conductivity of a metal. The 
latter gives a correct order of magnitude for v\\ and 

Having in mind all these facts one can readily 
write the equations of motion of the ion under a 
joint action of the AC electric field Eocosuit di- 
rected along z (normally to the surface) and a per- 
manent magnetic field directed parallel to the sur- 
face along y-axis, H y = H . 

1 dV 
M'z = -Mv{z)z + -QxH — - — h + QE (t) 
c oz 

(10) 

Mx = -Mv n (z)x + Q^ + C c (t) (11) 
c 

Here c is a light velocity. We included also into 
these equations the random forces £ x (z)(t) which 
reproduce the thermal fluctuations of the ion ve- 
locity. As usual, they have the white noise corre- 
lations of the following type 

< £*(ti)&(*a) >= 2u\\{z)k B T6(t 1 - ta) 

< Uh)Uh) >= 2v ± (z)k B T5(t 1 - t 2 ) (12) 

where ks is a Boltzmann constant and T is a tem- 
perature of the system. 
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In spite of the comparative complexity of these 
equations they can be analyzed analytically in the 
limit of the very small electric and magnetic fields. 
At small amplitudes E , H of external fields a dis- 
placement of the ion Sz = z — z e from its equilib- 
rium position is quite small and the first equation 
(10) could be chosen in a linear form 



Mz = -2Mzv -k e Sz + QE {t) 



(13) 



where k e is a force constant of a vibration of ion 
along the z-axis 

2 V 



k e — 



dz 2 



and 



u = 2v±(z e ) 



(14) 



(15) 



We neglected in (10) the thermal fluctuations (T = 
0) and the Lorentz force term. This is correct if we 
are interested in a linear contribution with respect 
to magnetic field into the velocity of the transport. 
However we have to keep this term in (11). 

The nonlinearity of Eq.(ll) is connected with 
damping term, which in accepted approximation 
could be written as 



- Mv (l - jSz)x 



(16) 



Here 7 is a derivative of the damping function 
v taken at equilibrium distance z e . Then, the 
Eq.(ll) took the following form 

Mx = -Mu (l - j5z)x + —z (17) 

c 

The periodic solution of (13) is 

Here lu c is an eigen-frequency of the physisorbed 
ion, u 2 = k e /M. Then the periodic solution of the 
second linear equation (17) could be written as the 
following 



x = Jo°° di[QH/cM] exp[-& ]6z(t - £) : 
exp {-7^0 ft_ £ Sz(t)cIt} 



(19) 



We are interested in the transverse velocity x aver- 
aged over the period of the vibration in the exter- 
nal electric field Eo(t). Making use the necessary 
averaging we have a final expression 

Q 3 E 2 H 7 lo 2 



v{uj) = 



1 



2cM 3 oj 2 + v 2 {uj 2 e - UJ 2 ) 2 + Auj 2 v 2 



(20) 



The function v(u>) has a resonance character and 
reaches its maximal value at w = lu c . Since 
vq <~ M _1 and lo 2 ~ M^ 1 the v(uj e ) depends quite 
weakly on the mass of ion. As a rule, v a <C oj c 
(damping is comparatively small) and 



v(0J e ) = 



Q 3 E 2 Hj 
2cM^Auj 2 e v 2 



(21) 



The behavior of v(u) as function of uo is given in 
Fig.l for three different values of mass. On the 
same figure we gave the results of numerical solu- 
tions of Eqs.(10,ll). 

Let us estimate the possible velocities v at rea- 
sonable experimental parameters 



M ~ 10 -g, 



24 f 

lO^sec" 1 , 



7 <~ 10 8 cm 1 , 
vo <~ 0.01w e 



and Q is an electron charge. Then 

tj(w e ) ~ W~ 3 E^H, 

where Eq and H have to be taken in emu sys- 
tem. For example, if H ~ 10 3 Gauss and Eq ~ 
100V/ cm, v(oj e ) ~ 1 cm/sec. 

In order to clarify the physical source of the di- 
rected transport of the particle we give a qualita- 
tive analysis of the equations(13) and (17). First of 
all we neglect the dependence of the damping coef- 
ficient v(z) on z, 7 = 0. Then, system of equations 
(13) and (17) became a linear one. The solution of 
this system can be written as following 

Sz = ao(u>) cos (cot + ip z (u))) 
x(t) = a (uj)b (uj) cos (tot + ip z (u>) + ip x {u>)) (22) 
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where a (oj) and <p z (uj) are modulus and argument 
of a complex function 

E Q 1 

M UJ 2 - LU 2 + l2v LO 

and bo(u>) and ip x (ui) are modulus and argument 
of a complex function . 

l/n + iuj 

A closed trajectory of the particle in the plane 
(5z,x) is an ellipse obeyed the following equation 
(Fig.2) 

(x - Szb cos(p x ) 2 5z 2 _ 
0565 sin ip x a Q 

If either ojl/^o ° r w/fo trends to zero the ellipse 
turns into a line. An angle i/> between axes z and a 
largest axis of the ellipse is proportional to bo(iv). 
Points A and A' on the Fig.2 are the points of re- 
turn of the trajectory in the re-direction. At points 
C and C the trajectiry cross the line 5z = 0, and 
dashed lines AB and A'B' are parallel to this line. 
In fact, the ellipse (23) is an attractor of the system 
(13),(17). 

Consider a small deformation of the trajectory 
taking into account the dependence of rj(z) on z. 
According to (16)the damping coefficient i](z) is 
higher at Sz < (7 > 0) and lower at 5z > 0. 
After switching on this dependence the trajectory 
of the particle between points C and B will be 
shifted a little bit to the right side. But this shift 
will be compensated completely by the shift to the 
left side on AC part of the trajectory. The same 
is true for AC and C'B' parts of the trajectory. 
Non-zero shift originates from AB and A'B' parts. 
Indeed the velocity on the AB part is positive and 
the trajectory will be shifted to the right side if 
we diminish the damping coefficient of the parti- 
cle. The velocity on A'B' is negative and if we 
increase the damping the trajectory will be shifted 
also to the right side. All in all it leads to the shift 
of whole trajectory to the right which is equivalent 
to the transport of particle in this direction. The 
qualitative consideration given above does not de- 



pends on details of the attractor and holds for any 
non- linear dynamics of the system at. 

Averaging the Eq.(17) over period of time we 
arrive to the following equality for the shift in re- 
direction 

Sx = 7 J Szxdt (24) 

The integral in (24) is nothing but an area of the 
ellipse S and we have a very simple formula 

Sx = jS (25) 

Eq.(25) could be used directly for the calculation 
of v(u) and of course it leads to the same result as 
(20). 

Two additional remarks to the main expression 
(20). If the external AC electric field is a sum 
of two electric fields with two different amplitudes 
and frequencies u\ and w 2 the effect is additive 

v = Ui(wi) +^2(^2) 

where fi(2)(wi(2)) are average velocities (20) with 
£0 = £-1(2) an d oj = Wi(2), correspondingly. More- 
over, any non-thermal electric noise in system leads 
to a directed transport of charged particles. Cor- 
responding average velocity could be obtained by 
integration of (20) over a frequency distribution of 
the electric noise. 

Now, if we have two weakly coupled particles 
with opposite charges and different masses Mi and 
M 2 the average velocity of their mass centre can 
be defined as 

_. , MivT(uj) - M 2 v5{lu) 

v(oj) = — — - 

y ' M 1 +M 2 
leading to some neutral current in system. More 
realistic situation with strongly coupled charged 
particles will be considered in Section 4. 

II. THE NLH EFFECT ON METAL 
SURFACE. 

The similar phenomenon takes play inside a 
metal on the boundary separating two media hav- 
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ing different conductivities. In this section we con- 
sider the metal-insulator or metal-vacuum bound- 
ary as simplest example of such interface. Con- 
sider dynamics of electronic gas in surface layer of 
metal under an influence of external AC electric 
field Eo(t) — Eq cos cot directed perpendicularly 
to surface of the metal (along z axis) and in the 
magnetic field H directed along the surface, axis 
y. At this geometry a density of electronic gas, 
p(z, t), is a function of normal coordinate z and a 
time t. There are two velocity components of the 
electronic gas: the normal component v(z,t) and 
the tangential component u(z,t) both of which are 
functions of only z and t. The electronic gas vi- 
brates perpendicularly to the surface. Due to the 
nonlinearity of this vibrations a normal velocity 
v(z,t) and density of electron gas p(z,t) are peri- 
odic functions of time with frequency u> having all 
overtones: 

V = v\(z) cos (u>t + ipi ) + t>2 (z) cos (2u>t + If 2 ) ■ ■ ■ 

P = Po+Pi {z) cos (u)t + ip 3 )+p2{z) cos (2ut + (fi) . . . 

The above amplitudes and phases of harmonics as 
functions of z are a subject of calculations. The 
Lorentz forces induced by magnetic field lead to 
tangential vibration of electronic gas with ampli- 
tude which is proportional to permanent magnetic 
field. Corresponding tangential velocity u(z, t) has 
the same form as v(z,t) with its own amplitudes 
and phases. The density of tangential current 
is j(z,t) = u(z,t)p(z,t). The time averaging of 
j{z, t) over period gives the DC component of the 
current. 

Since the width of a surface layer and all other 
characteristic lengths greater than atomic dis- 
tances one can readily use hydrodynamics ap- 
proach (see HI?]]) for description of system. In pla- 
nar geometry all variables are the functions of only 
one coordinate z and a periodic functions of t with 
frequency lo. The main equations look as following 




-pv 



rap. 



j du du 

p [m +v Wz 



dv \ 
V d-z)=- 


dp dz 


— Q{z) 

m 


-E (t)p- 
m 


He 
pu - 

mc 


oz z 



(26) 



He 



e d 2 u , . 

pv pu + 77— (27) 

mc rap oz z 



dp d 

m = -d- z {vp) 



Q(z,t) = 2q(z) -g(oo) 



(28) 
(29) 



Here p is a pressure of electronic gas in metal, p is a 
mobility of electron in metal, c is the light velocity, 
m is effective mass of electrons, rj is a viscosity 
coefficient of an electron fluid. A quantity 

q(z,t)=e (p(z',t) - p )dz' (30) 
Jo 

is a surface density of charge (together with pos- 
itive charge of a lattice) accumulated in a layer 
between z = and z. The boundary conditions 
for this system are quite simple. Both velocities 
v, u are equal to zero at the edge of metal (z = 0). 

We are looking for a periodic solution of the 
equations (22-25). The surface electric current di- 
rected along the surface and averaged over the pe- 
riod of time is 



ep(z, t)u(z, t)dz 



(31) 



The bar in Eq.(31) means the time averaging. 
The equation (28) can be rewritten via the q(z, t) 

as 



dq 
dt 



-ev(z, t)p(z, t) 



(32) 



From this equation we see immediately that the 
time averaged normal component of the current 

Then, using the Eq.(27) we can get the following 
formula for time averaged tangential component of 
the current 
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d 



d 2 u 



u»j x (z) = -e—[p(z,t)u(z,t)v(z)} + e V-g^ (33) 

The full surface electric current Ijj can be obtained 
by integration over z and looks as following 

epo e di] f°° dpdu 



-dz (34) 



u/,, mw M dp J dz dz 

where «oo(£) and Woo(i) are the normal and tan- 
gential velocities for large z. The constant 

drj 

K = m/ dp- 

has a dimension of a kinematic viscosity cm} j sec. 
At small electric and magnetic fields first non- 
vanishing contribution to the Ih can be obtained 
if we know q(z,t), u(z,t) and v{z 1 t) in a linear 
approximation. 

Taking the limit z — > oo we come to three 
equations for asymptotic values q x = q(oo,t), 

Uoc = u(00,t), Voo = v(00,t) 

d 



dt 



-pov l 



(35) 



j t v x = + ^E (t) - 20q % (t) (36) 



dt 



«oo = V l U) L - UiU 



(37) 



Here 



= e/mp, P = 2ire 2 /m and u>l = eH/mc 
is a Larmor frequency. Periodic solutions of these 
linear equations give the boundary conditions at 
z — > oo. And q(0, t) = by definition. Since one 
needs only the periodic solution of (26-29) the ini- 
tial conditions arc arbitrary. The system of equa- 
tions (26-29) has been analyzed both by perturba- 
tion theory with respect to external electric field 
and numerically using standard NAG routines. 

The periodic solution of linear equations (35-37) 
can be written in following form 



««,(*) = c e i?e[a(o;)e iw *] 

?«>(*) - [e c e /Lo^p Re[b{Lo)e luJt } 

jfeo(i) = e hc e Re[c{uj)e wt ] 



(38) 



Here we introduced a dimensionless functions 



a(u>) 



b(w) = 



(u) 2 -UJ 2 )+ iLUUJ^ 



c(lu) = a(u)) 



(uj 2 -uj 2 ) + iujui^ 



(39) 



UJ U + 



where eo and h are the dimensionless electric and 
magnetic fields 

eE eH 2 2ire 2 p 
e = , h = , uj a = 



m 



Using (26-29) as boundary conditions at an infi- 
nite z we arrive to the following solutions of the 
Eqs. (26-29) 

v(z,t) = c e Re[a{uj)e iut {l - e~ ik )] 

q(z, t) = p Re[b(uj)e iut (l - e^ k )} (40) 

u{z,t) = c e Re[c{w)e ibJt (l - e~ ik )] 

Here £ is a dimensionless coordinate, z = c^jui^ 
, and complex k = ki + ik 2 is defined from the 
equation 

(fci + ik 2 f = [2uj 2 s - uj 2 + iwu^/ul, k 2 > 0, 

(41) 

We are choosing a root which has k 2 > 0. Omitting 
intermediate calculations we give a final expression 
for Ih in following form 

/h = -^W7iM + ^72M; (42) 

7! and 72 are dimensionless functions of a fre- 
quency u> of the external field 

71 (u>) = a*(uj)c(u>) + c*(u>)a(u>); 

72M = [k 2 k*c*(Lu)b(uj) + h.c][k + k*}- 1 

Both of them have resonance character and reach 
their peaks at uj — lu s as it shown on Fig. 3. Max- 
imal value of them are 71 (u) s ) — 2u> ju) s and 
72(^5) = 1. For metals at low temperature as a 
rule uj 2 /oj 2 -C 1. At very high frequency lj ^> oj s 
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both 7i(w) and 72 (oj) decrease as l/ui 4 . At low 
frequency u> <C w M they behave as o> 2 : 71 (w) ~ 
vlv 2 /ui and 72 H = uP/u^. 

Thus, at low temperature main contribution to 
the surface current Ir comes from the second term 
of (38), depending on the viscosity of electronic 
gas. According to || the viscosity of electronic 
Fermi liquid in metal at low temperature could be 
estimated as 

T] = mp c e d(e F /T) 2 (43) 

where d is a lattice constant and £j? is a Fermi 
energy of a metal. 

Thus, at low temperature the surface current Ijj 
can be expressed as 

^ _ u 2 ep dc e f e F \ ( (iE\ ( pM\ 
4w 2 yk B T J y Ce J y c J 

It is easy to see that the temperature dependence 
of the effect is determined by the factor a 3 (T)/T 2 
(ct(T) is a conductivity of the metal). 

To demonstrate the value of the effect let us take 
material parameters of pure copper at low temper- 
ature, say, T=30 K: 

d = 1CT 8 era, c e = 10 s — , po = 10 22 cm" 3 . 

sec 

The conductivity behave as T~ 5 : 

a(T) = a R (T R /Tf 

where a R = 10 6 am is a room temperature con- 
ductivity and T R = 300 K. 

Finally, for the pure cooper at T = 300 K we 
have 

I H = -0.lv 2 ElH 

where v is expressed in MHz, Eq in a V/cm, H in a 
Tesla, Ih is in a \iAjcm. At Eq = 10 V/cm, H = 
0.01 T, v = 10 MHz one has I H =~ 1.0 pA/cm. 



III. NLH CURRENT IN AN INJECTION 
LAYER OF A SEMICONDUCTOR. 

In the previous Section we have considered the 
dynamics of electron gas inside a metal on the 
boundary metal-semiconductor. If a work func- 
tion for transition of electrons from metal side to 
the semiconductor is not too high electrons of the 
metal can be injected into the semiconductor creat- 
ing an injection layer of electrons inside of the semi- 
conductor. The injection layer can give a contri- 
bution to the NLH current on the interface metal- 
semiconductor. An estimation of the NLH current 
in this injection layer is the subject of present sec- 
tion. First of all, let us consider the equilibrium 
electron density in the semiconductor side. Free 
energy of electron gas per unit area of the surface 
has the following form || Liebsch] 

oc 00 oc 

F = J f(p,T)dz + 27re 2 / p(z)dz J zp(z)dz- 
00 
00 00 00 
e 2 7r J J p(z)\z ~ z'\p(z')dzdz' + ip J p(z)dz (44) 
00 

where / = (p(z),T) is a density of the free en- 
ergy of injected electron gas, which is a function of 
electron density p(z) and temperature T. At small 
density p(z) the / = (p, T) is a free energy of ideal 
classical gas. 

f(p,T) = Tp\np/pT 

We use the units where Boltzmann constant 

fee = 1- ^0 is.a work function for transfer electron 

from metal to the semiconductor. Minimization 

of F with respect to p{z) leads to the well known 

equation for equilibrium density p e (z) 

ft f 00 00 

+ 2-Ke 2 z J p e {z)dz - e 2 2iT J p e (z')\z - z'\dz' 

VP 

oc 

+ 2vre 2 J zp e (z)dz + f = 


(45) 

For Boltzmann electron gas the solution of (45) 
looks as following 
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Pe(z) 



Zdqo 



Here = p{z)dz and 

2^e 2 



I) 2 



(46) 



Zd 



T 



is full injected charge per unit of area of surface. 
An equality of chemical potential of the metal and 
the semiconductor on the surface of the metal 



2i 
dp 



P=Pe(0) 



gives 



and 



p e (0) = p T e 



-*/T 



(PT/z d 



,l/2 e -*/2T 



(47) 



(48) 



(49) 



In fact, the equations (45-49) summarise briefly the 
simplest theory of Schottky barrier on the metal- 
semiconductor interface ]9| . Under action of time- 
periodic external electric field Eq sin cot the non- 
linear oscillations of the density p(z,t) takes place. 
A system of equations describing these oscillations 
looks very much like that of the metal-vacuum in- 
terface equations (26-29) , but the background pos- 
itive charge po in Eq.(29) has to be taken zero. 

Another important difference between the 
metal-insulator and metal-semiconductor cases is 
connected with boundary conditions for velocities 
v(z,t) and u(z,t) at z = 0. There are two limit 
cases. In the first one, the exchange of charges 
over the surface is quite slow which means that to- 
tal charge in the semiconductor is conserving un- 
der action of the AC electric field and v(0,t) = 0, 
u(0, t) = at the boundary. If we neglect the vis- 
cosity a total NLH current vanishes. 

Let us outline briefly the solution of the hydro- 
dynamics equations (26-29) in this case. At small 
electric Eq and magnetic H fields and small vis- 
cosity of electron liquid the first equation (26) can 
be rewritten as 



where 
and 



,dp 

'dz 



+ pMp(z)-p 00 ) + eE(t)p (50) 



p(z,t) = / p(z)dz 
Jo 



"dp 2 



(51) 



(52) 



For Boltzmann gas c 2 e = y^T/m . The equation of 
continuity can be rewritten as 



dq 

di = - VP 



(53) 



In the equation for u(z,t) we have to keep small 
terms connected with the magnetic field (Lorentz 
force) and viscosity. Then, this equation can be 
rewritten in the following form 



dup d 



— — + -^-vup = uj L vp - ujfiup H (54) 

at oz m dz 

Having in mind that we are interested in the full 
surface current averaged over period of oscillation 
of external electric field 



Ih= e / pudz 



(55) 



(56) 



we arrive to the expression 

e f°° 77 d 2 u , 
Ih= — / — — ^dz 

w M Jo m dz 

Here again the bar means the averaging over period 
of oscillation. We took into account also that 



jr. 



vp= 



(57) 



due to the second equation. 

As a rule, the viscosity of electronic gas 77 de- 
pends on the density of electron and, thus, depend 
on z . It is more convenient to rewrite (56) as 



Ih= — 



oj,. m 



^fj. •<" jo 

Further we shall show that 

du 
dz 

z=0 



,di] dp du ^ 
dp dz dz 



(58) 
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At small viscosity and small uj/uj^ (low frequency 
of external field) 



U = V 



and we can express the In via normal velocity 
v(z, t) and density p(z, t) found from equation (50) 

euj L f°° drj dpdu 

Ih= 2 / (7T^7Ta~ dz 59 ) 

moj^ J op Oz oz 

The viscosity of electron gas is unknown function 
of density of the gas. Father, for estimation, we as- 
sume that is — ^ is approximately equal to kine- 
m op 

matic viscosity of electronic gas, 

1 Or/ 

K= — JT> 
m op 

taken at equilibrium surface density of electron 
gas, p eg (0). Thus (60) turns into the following ex- 
pression 

Ih= ~—^L **TzTz (60) 

which we will use for estimation of nonlinear Hall 
surface current. 

The solution of the eqs. (50) will be performed 
by theory of perturbation with respect to a small 
amplitude of electric field E . After introduc- 
tion of dimensionless time r = cut , coordinate 
'•• i" , . . .. , , . ' l-'u 



function X(z, t). 



, electric field eo 



mCgy/UJUJ^ 



and new 



Q(z,t) = <?oo(l - aX(z,t)) (61) 
the (50) turns into the following equation 



dX OX 



!Tr=W + 271 dZ +6 dZ COST (62) 
which has to be solved with boundary conditions 

A(£) -> at £ -> oo 



A(0) = 



a 



Here we use the following notation 



\M^c e m Ce 
a = \ —>P = —F= 



dX dX 

Note that we are looking for the periodic solution 
of (63). If we make a substitution 



a =-<!>/, 

the nonlinear eq. (63) turns into linear one 

dtp d 2 ip d<p 
— — = — tt + e cos t— - 

with the boundary conditions 

<p(£) — > oo at £ — > oo 



(63) 



(64) 



d<p 



-p(0) 
a 



The perturbative analysis of linear Eq. (65) is 
rather trivial and we give only the final result 



X(£,t) = XoiO+eX^t) 



A = 



Ai = e 



1 da(0 



where a(£) has the form 



o(0 



a(0 = i 



1 



1 + fca 



-fc? 



- 1 



t + 1 
V2 



(65) 



(66) 



(67) 



(68) 



Omitting intermediate calculations we give the ex- 
pression for Ih in the following form 



bui L goo a 



Z 3 7o (£ + a) 



/°°77^4 V ^)^ ( 69 ) 
Jo (£ + a) 



where 



•<«-(^t)/(t: 

since is proportional to v 7 ^ ,the function i?(a) , 
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R(a) = 



1 



(70) 



is actually a function of frequency of external 
electric field ^Numerical calculation shows that 
R(0) — —0.21. Collecting all multipliers we get 
the following expression for the Ih 

/ \ 2 

(71) 



As an example let us choose reasonable parameters 

i9 1 q cm 2 
Poo = 10 12 t, k = 10 3 , 



/i = IO 5 , v = IO 7 — . 

V sec sec 

Then, Ih can be written as 

I H = lO~ 5 HuE$ 
llA 

where Ih is expressed in a , H in a Oe, in 

cm 

a V/ cm and w in a MHz. 

IV. MAGNETOPHORESIS OF NEUTRAL 
PARTICLES ON THE SURFACE. 

Consider now the neutral particle (atom, 
molecule or any nano-size particle) located on the 
metal surface under joint action of normal AC- 
electric and tangent permanent magnetic fields. As 
it has been shown in Section l the driving force 
which causes the directed transport of particle on 
the surface depends on a mass, charge and details 
of an interaction between particles and the surface 
of metal. If negative and positive charges coupled 
in a neutral pair as it takes place in atoms and 
molecules then the driving forces acting on positive 
and negative charges are directed on opposite sides, 
but their values, generally speaking, can be differ- 
ent. Thus, the resulting force acting on the par- 
ticle is not necessary zero leading to the directed 
transport of the particles. As it is well-known the 



Lorentz forces acting on the neutral particle in a 
magnetic field H can be written as 



F, = i 

c 



<9D 

~dt 



x H 



(72) 



where D is a dipole moment of the system of the 
charged particles. On the other hand, the dipole 
moment D of the neutral particle in an external 
electric field Eo cos(wi) is an algebraic sum of a 
permanent dipole moment of the particle Do and 
the dipole moment Di(i) induced by an external 
time-periodic electric field Eocos(ui). Thus, we 
have 



D(t) = D + a(w)E cos(u^) 



(73) 



where a(u>) is a frequency depending polarisibil- 
ity of the particle. As a rule, for neutral particle 
absorbed on a metal surface vector Do is directed 
normally to the surface. The potential of interac- 
tion between the dipole particle and metallic sur- 
face V(z), consists of two terms ||. 
The Coulomb attraction potential 



V id {z) 



D 2 

4z 3 e 



(74) 



and a repulsive exchange potential V e d{z), which 
we again (see Section l) choose in an exponential 
form 



V ed {z) = c d a d exp 



Zd 



a,] 



(75) 



Here a d is a radius of exchange forces on the 
metallic surface and z d is an equilibrium distance 
between the particle and the surface, c d is deter- 
mined from an equilibrium condition 



V{z) = V id {z) + V ed (z), 



dV_ 

dz 



(76) 



As in the case of ion-metal system the friction 
coefficient, i/, for the dipole particle located on a 
metal surface depends on the distance z (see ||). 
At large distances z, the v(z) is proportional to 
l/z 3 . 
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Similar to Section 1 we choose this dependence 
in a form 



v{z) = v [l -j d (z - z d )\ 
dv{z) 



(77) 



where Vq = v{z d ) and 7^ = 



dz 



At small deviation from the equilibrium position 
of the particle the equations of motion for its grav- 
ity center can be written as following 

Mi = -Mv(z)z - u) 2 d M(z - z d ) - Q d E (t) 
H 3D 

Mi = -Mv{z)x H — (78) 

Here uo d is an eigenfrequency of particle vibra- 
tion on the surface and 

3Doo;(w) 



Using (73) we get 



dD 

= —u)a(uj)E sin(wt) 



(79) 



The system of equations (78) is complete analog 
of that in the Section 1. Omitting the details of 
the calculation we give the expression for tangent 
velocity averaged over a period of oscillation of 
the external effective electric field 

± = V (w) = \ [sbMriM + ffiMsSM] (so) 

where complex functions go(u)) and g\(u)) have 
form 



ffoM 



QdEo 



M{u 2 d - uo 2 + iv u) 
E a{uj)Hiuo 



Mc(iuj + v ) 
Final expression of the v(u) has a form 
3L>o a 2 (uj)E 2 HLU 2 



(81) 
(82) 



v(uj) = 7 



zjM 2 c 



+ vt (lu 2 - uj 2 Y + V 2 UJ 2 
(83) 



At small u> the velocity v(w) is decreasing as oj 2 
in complete analogy with corresponding velocity of 
charged particle. However, in difference with the 



Section 1 v(ui) changes its sign at u> = ^/u> d + Vq. 
The general form of v(lu) as a function of frequency 
is presented on Fig. 4. 

If the dipolc moment of the particle D is pro- 
portional to its volume (as it takes place for ferro- 
electric medium) then v(u>) does not change much 
with the increase of the radius of the particle and 
its mass. 

If the external electric field Eo(t) is a sum of two 
harmonics with different frequencies uj\ and iOi the 
resulting velocity is additive, 



v{uS) = wi(wi) + v 1 (uj 2 ) 



(84) 



For atomic size particles (atom or molecule) the 
maximal velocity could be estimated as 



v = 1Q- 12 E 2 H 



(85) 



where v is given in cm/sec, E in V/cm and H in 
Oe. 



V. NLH EFFECT PHENOMENOLOGY. 

In previous Sections we have shown that the 
joint action of the permanent magnetic and AC- 
electric fields leads to the direct transport of par- 
ticles located nearly surface. An important role in 
consideration plays gradient of the damping con- 
stant the system metal-particle. Generalizing all 
examples of the previous Section we can assume 
that the NHL effect takes place in a bulk of the 
material if there exists a non-zero gradient of con- 
ductivity tensor ji could be chosen as 



ji = \{u)—-^-E k Eieji m H„ 

oz k 



(86) 



where E k is a time derivative of an external electric 
field and H is magnetic field. In a case of isotropic 
medium 07, = S^a the density of current is 



j = A(VctE)[E x H] 



(87) 



Here A is a phenomenological coefficient depending 
on frequency and material parameters. 
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To demonstrate the physical consequences of 
consider the cylindrical system shown on 
Fig. 5. The time periodic voltage is applied be- 
tween inner and external surface of the cylindrical 
sample. The permanent magnetic fields H is di- 
rected along main axes of the cylinder. We assume 

also that the conductivity of the metal decreases 

dcr 

along the radius of the sample, — < 0. 

or 

In the chosen cylindrical geometry there exist 
only angular component of the current density 
3<p ( r ) which according to (88) can be expressed as 



Jv (r) = AM^E H 
Or 

Now, let us assume also, that the conductivity a 
, u> and thickness of the sample arc chosen such a 
way that the electric field E(t) penetrates over the 
sample. Then the additional magnetic field H a d(r) 
induced by this current can be written as 

47T [ r 47T 

Hi(r) = / Jv (r)dr = — (r - (89) 

c J ri c 

where n is an radius of the sample. There ex- 
ist also diamagnetic contribution to the induced 
magnetic field which is proportional to diamag- 
netic susceptibility, x >°f the material 



H d = X H 



(90) 



Summing up both terms we come to the conclu- 
sion that the full induced magnetic field is positive 
if 



- A— E 2 — l d > \ X \ 

or c 



(91) 



Actually, the condition (91) is a condition of sta- 
bility of the system to the generation of magnetic 
field in the sample under action of time periodic 
electric field E(t). The details of phenomenon will 
be considered in the future. 

The main obstacle for observation of the NLH ef- 
fect in AC electric field seems to be connected with 
an absorption of ultrahigh frequency electromag- 
netic field into material. The time-periodic com- 
ponents of Hall current which is linear with respect 



to external electric field is greater than a perma- 
nent component of the NLH current and may lead 
to heating up of the sample. On the other hand, at 
low (or ultralow) temperature the absorption and 
warming up could be made small enough for such 
observation. 
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Captions to Figures 

Fig. 1. Frequency dependence of the velocity 
(20) of a single charged particle at vq = 0.1w e 
for three different masses m = M/M = 0.25, 
M/M a = 1.0 and M/M = 4.0. 

Fig. 2. An attractor of the system (13-17). The 
A and A' are points of return in x direction. Ar- 
rows on the trajectory show the direction of the 
particle rotation, is an angle of inclination of 
the ellipse. 

Fig. 3. Frequency dependence of 71 (w) and 
72 (w) at cj m = 0.1w s . 

Fig. 4. Frequency dependence of the velocity 
(83) of neutral particle at vq = O.lo^- 

Fig. 5. The principle set up for observation of 
NHL effect. Electric field E (t) is directed along 
the radius of the cylinder sample. 
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